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Abstract 

Wave functions of bounded quantum systems with time- independent 
potentials, being almost periodic functions, cannot have time asymp- 
totics as in classical chaos. However, bounded quantum systems with 
time-dependent interactions, as used in quantum control, may have 
continuous spectrum and the rate of growth of observables is an issue 
of both theoretical and practical concern. 

Rates of growth in quantum mechanics are discussed by construct- 
ing quantities with the same physical meaning as those involved in the 
classical Lyapunov exponent. A generalized notion of quantum sensi- 
tive dependence is introduced and the mathematical structure of the 
operator matrix elements that correspond to different types of growth 
is characterized. 

Keywords: Quantum Lyapunov exponent, quantum complexity, quantum 
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1 Introduction 

Bounded classical systems that are chaotic, display exponential growth of 
initial perturbations and other interesting long-time asymptotics, like expo- 
nential decay of correlations. In contrast, quantum Hamiltonians of bounded 
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systems with time-independent potentials, having discrete spectrum, their 
wave functions are almost periodic functions. For this reason the work on 
"quantum chaos" has shifted from consideration of long-time properties to 
the statistics of energy levels of quantum systems with a chaotic classical 
counterpart (for a review of recent work see [||] and references therein). 

However, quantum systems with bounded configuration space but time- 
dependent interactions (for example particles in an accelerator subjected to 
electromagnetic kicks or the systems used in quantum control) may have 
continuous spectrum. Therefore the estimation and control of the rate of 
growth, of the perturbed matrix elements of observables, becomes an issue 
of both theoretical and practical concern. 

In classical mechanics the most important asymptotic indicator of chaotic 
behavior is the Lyapunov exponent (an ergodic invariant). Therefore a nat- 
ural first step to discuss rates of growth in quantum mechanics seems to be 
the construction of a quantum Lyapunov exponent. After a few heuristic 
attempts by several authors (see references in ||) a satisfactory construction 
has been achieved in the sense that the phase-space observables that are 
used are exactly the same in classical and quantum mechanics. It uses the to- 
mographic formulation which describes conventional quantum mechanics by 
a set of marginal probability densities || . Then, the only difference between 
the classical and the quantum exponent lies in the time evolution dynamical 
equation. In Section 2 we recall this result. Translating it from tomographic 
densities to traces of operators, it turns out that the quantum Lyapunov ex- 
ponent measures the rate of growth of the trace of position and momentum 
observables starting from a singular initial density matrix. 

A positive Lyapunov exponent corresponds to exponential growth of these 
traces. However, the same quantities may serve to characterize other types 
of growth, leading to a generalized notion of quantum sensitive dependence. 

There are examples where exponential rates of growth (as in classical 
chaos) are also found in quantum systems (Sect. 3). However, in many other 
cases, quantum mechanics seems to have a definite taming effect on classical 
chaos. Therefore, a generalized notion of quantum sensitive dependence, 
corresponding to rates of growth milder than exponential, might be of interest 
to classify different types of quantum complexity or to characterize the degree 
of accuracy achievable in quantum control. 

Quantum sensitive dependence is discussed in Sect. 4, as well as the math- 
ematical structure of the operator matrix elements, in the spectral represen- 
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tation of the trace, that corresponds to each type of growth. A convenient 
unified framework to discuss these matters is the space of ultradistributions 
of compact support and their Fourier images [|J 0. 



2 The quantum Lyapunov exponent 

As a first step we will rewrite the results of Ref. @] using operator traces. 
In Ref. the quantum Lyapunov exponent along the phase-space vector 
v = (viv 2 ) is shown to be 



A„ = lim — log 



( 



) 



where M t (X, /i, u) 

M t (X, n,v) = j IT (X, /i, v, X', //, u', t, 0) M (X, ^ v) dX' n d^ n dv' n (2) 

is the time evolved tomographic density ||, starting from the initial condition 
M (X', //, !/') = ((«! © fi' + v 2 © v') . V x >) S n {X' - fi'qa - u'p ) (3) 
a © 6 being defined as 

(a © 6) f = aA 

When II is the classical propagator, A„, as defined in eq. ([!]), coincides with 
the usual classical Lyapunov exponent constructed from the tangent map. 
The only difference between classical and quantum Lyapunov exponents lies 
in the dynamical law of the propagator, thus insuring that we are dealing 
with quantities with the same physical meaning. 

For a system with Hamiltonian 

H = £ + V{q), (4) 

the evolution equation for the quantum propagator of the tomographic den- 
sities is 

§ - n • VJI - V X V (q)m(u® VxIT) 

+1 E^-ir 1 (f) 2n+1 Vil ^r @ (" © Vx) n ■ ■ > ® v x y t2n+i n (5) 
= o 
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with initial condition 



lim IT (X, ji, v, X', p', v', t, t ) = 5 n (X- X') 5 n (p - p') 5 n [y - u') (6) 

t— »t(j 



reducing for h = to the classical evolution equation. 

In the tomographic formulation, classical and quantum mechanics are 
both described by a set of positive probability distributions M t (X, p, u), 
the H— deformation appearing only in the time-evolution. It is this fact that 
allows the notion of Lyapunov exponent to be carried over without ambiguity 
from classical to quantum mechanics. However, to relate the Lyapunov expo- 
nent to the behavior of operator matrix elements and the spectral properties 
of the Hamiltonian, it is more convenient to rewrite it as a functional of the 
density matrix p (x, x ) . The first step is to consider the Fourier transform 
G t {p, p) of the tomographic density M t (X, p, v) 



G t (p, v) = G t (1, fi, v) = J d n X e tXml M t (X, p, v) 
and perform the integrals in (|]) to obtain 
A 



>'2 



lim - log 

t^oo t 



V^G t (p, v) 
V„G t (p, v) 



|/i=i^=0 
fi=v=0 



(7) 



(8) 



Now, using the relation between the tomographic densities and the density 
matrix [0, namely 



G t {p, v) 



1 

2^ 



one easily obtains 



A 



>'2 



d n X d n pd n xd n x' ei Xml - p '( x - x ')) Pt (x,x) (9) 
5 n [X - p® ( X X ) +v ®p 



(10) 



lim - log 

t^oo t 



Tr {p t x} 
Tr {ptp} 



the density matrix at time zero (corresponding to Mo (X ,/x',z/') in eq.@) 
being 
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Po[x,x 



—e 



ipo* I x—x 



[vi • V) 5 n [ q 



X + X 



+ iv 2 



go 



X + X 



(11) 

Eq. (|lO|) means that the quantum Lyapunov exponent measures the expo- 
nential rate of growth of the expectation values of position and momentum, 
starting from the initial singular perturbation p . This is a rather appealing 
and intuitive form for the Lyapunov exponent. That the quantum Lyapunov 
exponent should have a form of this type had already been proposed in 
Ref. 0, based on qualitative physical considerations. What is not obvious, 
though, without the tomographic formulation, is that this is the form that 
corresponds exactly to the same physical quantity as the classical Lyapunov 
exponent. Also non-obvious is the specific form that the initial singular per- 
turbation p should take. 

Using the time-dependent operators in the Heisenberg picture 

XH (t) = WxU 
Ph (t) = WpU 

one has an equivalent form for A-> 



(12) 



A 



Vl 
V-2 



lim - lo£ 



Tr' {p x H (t)} 
Tr' {poPH (t)} 



(13) 



where we have also defined 

Tr {p x H (t)} = Tr {p x H (t)} /Tr {p x H (0)} 

Whenever PqXh it) is a trace class operator, the term corresponding to Tr{p XH (0)} 
has no contribution in the t — > oo limit. On the other hand, by taking the 
appropriate cut-off and a limiting procedure, the above expression may also 
make mathematical sense even in some non-trace class cases. 



3 An example: Kicked motions in the torus 

Let x\, X2 G [—iv, 7i") be coordinates in the 2-torus T 2 with conjugate momenta 
Pi,P2 and the dynamics be defined by the Hamiltonian 

H = H + J2 V (x,p)8(t-nr) (14) 
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x G T 2 , p G R 2 and, in particular, H = ^ . Let 

U= C 2 ([-7T,7T),C/ 2 X) (15) 

Physical observables should be self-adjoint operators. Therefore the domain 
D (xi) of Xi is 

D(x, i ) = {fEH} (16) 

and the domain of Pi 

D (pO ={/GW|/(i { = -7r) = / {Xi = vr)} (17) 
A convenient basis of vectors in D (p^ is 

H = {(x\q) = -^=e^ x \keZ 2 } (18) 
The Floquet operator associated to the periodic Hamiltonian H is 

Up = U U K (19) 

with Uq = exp (iH t) and 

U K = exp(iV(x,p)) (20) 

We will consider different types of kick potentials. Whenever V(x,p) is a 
function of x or p alone, any differentiable function will generate an unitary 
operator Uk operating in the basis flTE|). However for kicks of the electro- 
magnetic type, V (x,p) = | [x-iPi + PiXi), because 



D [ \ (xiPi + PiXi) ) = D (xiPi) n D (piX-i) (21) 



2 



{fen\f(x l = -n) = -f(x t = n)} 



\ {xiPi + PiXi) does not generate a continuous unitary group in (|l"8"l) and only 
a discrete set of kicks will be acceptable, namely 

U K = e*p(±(x.A.p+p.A.x)) (22) 
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exp (A) being a 2x2 matrix with integers entries and determinant one, the 
last condition resulting from 



e ±(x9A»p+pmA*x) _ e x»A*p e ±TrA 



Using the momentum basis (|T8[) the initial density matrix p of eq.(|TT|) may 
be written 



( vi • + v 2 . £-) V |po - k) (po + k\ (23) 



Po = -2 

and in a position (generalized) eigenstate basis 

Po = -^( Vl .^ + v 2 *^-) [ d 2 x\q + x)e^(q -x\ (24) 
The corresponding traces, needed to compute the Lyapunov exponent, are 

TV { Po U } ( I ) U) = -2^. A +0a . |_) £ (p. + fc| ( ^ 

(25) 

and 



u t bo - fe> 



TrUoUH X p )u t } = -4.(, 1 .|- + , 2 .|-, (26) 



Another form, that will be used later on, is a spectral decomposition using 
the eigenmodes of the Floquet operator. For discrete spectrum 

Tr {pbflj (p)^*} = EW P° W ( p ) W c^-^* (27) 



;r) 



and in general 



Tr {po^/ (l) U ] = I dE»dE vPo (/i, „) ( ) c -*(^)* ( 28 ) 

Three types of potentials will be considered: 
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(i) V(x,p) = 

This is just free motion on the torus with U t = exp (i^-tj . Then 

and from @ it follows that Trjpo^ ( * J U t \ is a constant independent 

of time, implying = 0. A similar conclusion would be obtained analyz- 
ing the spectral decomposition because the spectrum of the Floquet being 
discrete in this case the right-hand-side of Eq.(^7p is an almost periodic func- 
tion. 

Free motion having an irrelevant effect on the computation of the Lya- 
punov exponent, we restrict ourselves, for simplicity, to the resonant case, 
r = 4nm, m G Z, in the next two examples. 

(ii) V (x,p) = ag(x) and r = 47rm, m G Z 
From 

(Po + *l Ul ( I ) U t bo - *> = <po + k\ ( I + t _ aVg {x) ) \p - k) 

it follows that Tr i^poLf} ^ X j U^j grows at most linearly with t, implying 
also A„ = 0. 

In this case the Floquet operator spectrum is continuous but the kernels 
x (is, p) and p (is, p) being 

x (v, p) ~ 5 (E v - Efj) 
p(is,p) ~ 



we obtain the same conclusion from the spectral representation (28) 

(iii) V (x,p) = \(x • A* p + p • Am x), r = A-nm, m G Z and M 
exp (A) = 



2 

1 1 ' 
1 2 

Let t = nr. From 



(p + k\ U\xU t \p -k) — (M~ lT y (po + k\x \p - k) 
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and 
with 



(p + k\ U}pU t \ Po -k) = M n (p + k\ p |p - fc) 



M r 



M~ 



U-Zn+l + w 2n-l _^-2n + ^2n 

^2n+l + u -2n-l _J2n + ^-2n 

-cj 2n + cj~ 2n u; 2 ™- 1 + uj- 2n+l 



and a; = ~ (l + Vo) , we conclude working out Eq. that in this case there 
is a non-zero Lyapunov exponent A^ = 2 logo;. 

It is instructive to find out how the same result may be obtained from the 
spectral representation fl2"8|). In the resonant case (r = 4nm) the eigenstates 
of the Floquet operator are @ 



1 oo 

\E a ) = -=Y" e ~ ian \M n P) 

V 27T — 

v n=— oo 

with eigenvalue exp (— ia). From 

(Ep\UtpU t \E a ) = M n (E^\p\E a ) 

= e-^ Ea - E » )nT (E^p\E a ) 

and the corresponding equation for x, it follows that the kernels x (u, p) and 
p (is, p) are linear combinations of 

5 (^E a -E^- 1 - log A fc 

Afc , k = 1, 2 being the eigenvalues of the matrix M. It is the complex shift 
in the argument of the delta that converts the complex exponentials in the 
spectral decomposition (p8|) into an exponential growing quantity. 

Both in case (ii) and (iii) the Floquet spectrum is absolutely continuous. 
Nevertheless the rate of growth of the traces is quite different. These exam- 
ples suggest that the critical role is actually played by the analytic nature 
of the phase-space operator kernels x (is, and p(z/, p). This will further 
clarified in Sect. 4. 
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4 Sensitive dependence in quantum mechan- 
ics 



Let us denote the dynamical variable appearing in Eq.flTBD as 



A(t) 



(29) 



Tr' {p x H (t)} 
Tr' {pqPh (t)} 

In the three examples studied in the preceding section this dynamical quan- 
tity shows no growth in the free motion case, polynomial growth for space- 
dependent kicks and exponential growth for the electromagnetic-like kicks. 
The second well as the study of the standard map in Ref. @ , clearly 

show the taming effect that quantum mechanics has on classical chaos. Nev- 
ertheless A (t), as defined in Eq . (|29f) , is the quantum observable that corre- 
sponds to the notion of separation of nearby trajectories in the classical case. 
Therefore, for example, a polynomial growth of this observable signals an 
higher degree of dynamical complexity that no growth at all. Therefore, in 
view of the widespread quantum suppression of exponential growth, it makes 
sense to characterize different degrees of quantum dynamical complexity by 
a more general notion of sensitive dependence. We define: 

Definition Quantum dynamics is sensitive- dependent in the support of 
p if for any T and M > 0, there is at>T such that (A (t) /A (0)) > M. 

The above definition allows for rates of growth slower than exponential 
and even for oscillations of the ratio A (t) / A (0). It only requires it to be 
unbounded. 

In the examples of the preceding section, free motion is not sensitive- 
dependent, whereas the x— dependent kicks in 3.(ii) are polynomial sensitive- 
dependent and the non-local (electromagnetic) kicks in 3.(iii) are exponential 
sensit i ve- dep endent . 

A precise mathematical characterization of when each type of sensitive- 
dependence is to be expected, is possible. This uses the well-known space of 
ultradistributions of compact support |§] || and the corresponding Fourier 
image. 

Let X n = {z : z G C, \z\ > n} and B n be the Banach space of complex 
functions analytic in X n and continuous in X n with norm 

= sup \<j>{z)\ (30) 

z€X n 
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The space of ultradistributions of compact support IA C is the inductive limit 
of the spaces B n . Its dual is the space of entire functions. An important 
subspace of U c is the space of distributions of compact support V c , the cor- 
respondence being established by the injective (but not surjective) mapping 
(the Stieltjes transform) 



Whenever / is a distribution of compact support, the ultradistribution Sf € 
U c vanishes at infinity. 

On the other hand the Fourier transform establishes a correspondence 
between the space of ultradistributions of compact support and the space of 
functions of exponential growth. 

An entire function \I/ (z) is said to be of exponential growth if and only if 
there are constants a and (3 such that \ty (z)\ < ae^' z ', \/z E C. The vector 
space of these functions will be denoted by 7i e . The Fourier transform 



with e IA C is a bijective linear map of IA C over Ti. e . 

On the other hand the restriction of T to the subspace T> c establishes an 
isomorphism between T> c and the subspace of 7i e consisting of entire functions 
with polynomially bounded growth on horizontal strips around the real axis. 

Now noticing that, in the energy differences (E^ — E v ) variable, the inte- 
gral in Eq. (|28|) is the Fourier transform of the kernels 



we conclude: 

Proposition A necessary condition for exponential sensitive dependence 
in quantum dynamics is that the kernel K, (/i, v) as a function of the energy 
differences — E v be a member ofU c /V c . If the kernel belongs to V c then 
there is, at most, polynomial growth. 

Cases (ii) and (iii) in Section 3 are examples where the kernels belong in 
the first case to T> c and in the second to U c /T> c . 




(31) 




(32) 




(33) 
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